Spectral statistics of a quantum interval-exchange map 
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Curious spectral properties of an ensemble of random unitary matrices appearing in the quanti- 
zation of a map p— >p-|-Q, g— >g-|- f{p + a) in Q are investigated. When a = m/n with integer 
co-prime m, n and matrix dimension A'' ^ oo is such that mN = ±1 mod n, local spectral statistics 
of this ensemble tends to the semi-Poisson distribution |3 with arbitrary integer or half-integer level 
repulsion at small distances: R2{s) — > s'' when s — » and /? = n — 1 or n/2 — 1 depending on 
time-reversal symmetry of the map. 
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Random matrices appear naturally in many different 
physical and mathematical problems ranging from nu- 
clear physics to number theory (see e.g. jS]). Though 
many different ensembles of random matrices were con- 
sidered, it is well accepted 9] that the level repulsion 
(i.e. the small-s behaviour of the two-point correlation 

function i?2(s) s^) is determined only from symme- 
try arguments: /3 = 1,2,4 for, respectively, real symmet- 
ric, complex hermitian, and self-dual quaternion matrices 
(plus, of course, /3 = for diagonal matrices with inde- 
pendent entries). Dyson's arguments leading to these 
values are so robust that only rarely one considers seri- 
ously the possibility that certain models may have differ- 
ent values of (3 (see e.g. ). 

The purpose of this Letter is to demonstrate that a 
quantization of an interval-exchange map leads to a ran- 
dom unitary matrix ensemble which, under certain con- 
ditions, reveals level repulsion with arbitrary integer and 
half-integer exponent /3. To the authors knowledge it is 
the first example which shows naturally level repulsion 
different from standard values (but see 

According to 0] a suitable quantization of the map 



f{p 



mod 1 



(1) 



(where a is a constant and f{q) is a certain function with 
period 1) leads in the momentum representation to the 
following N X N unitary matrix {k,p — 0,1, N — 1) 



Mkp = e} 



1 - e 



2-KiaN 



7V(1 - Q2Ki{k-p+aN)/N-^ 



(2) 



where $fc = -2nNF{k/N) and F'{p) = f{p). Its eigen- 
values have the form A(j) = e"''*--'-' with real called 
eigenphases and we are interested in their distribution. 

To define an ensemble of random matrices it is conve- 
nient to consider $fc in (|2Jl not as derived from a func- 
tion F{k) but as independent random variables More 
precisely, the following two cases are considered. In the 
first one (called non-symmetric matrices) all $/£ with 
k = 0, . . . N ~ 1 are independent random variables dis- 
tributed uniformly between and 27r. In the second case 



(called symmetric matrices) only a half of is indepen- 
dent. The others are obtained from symmetry relations 
^N-k = $fc- These two cases correspond to classical 
maps without and with time-reversal invariance. 

When a is a 'good' irrational number (e.g. a — \/5/2) 
numerical calculations show that spectral statistics of 
symmetric (resp. non-symmetric) matrices is well de- 
scribed for large N by the standard Gaussian orthogonal 
(resp. unitary) ensemble of random matrices as it follows 
from usual symmetry considerations. 

For rational values of a = m/n with integer co-prime 
m and n only a finite number of different momenta ap- 
pear under the iterations and classical map can be 
identified with an interval-exchange map. 

Our main result in this case is the following. When 
a = m/n and matrix dimension A'^ — s- cx) is such that 



mN = ±1 mod n 



(3) 



local spectral statistics of matrix ensemble (01 is de- 
scribed by the semi-Poisson statistics characterized by 
a parameter (3 related with a as follows 



/3 



n ~ 1 for non-symmetric matrices 
n/2 — 1 for symmetric matrices 



(4) 



The general semi-Poisson statistics was investigated in [Sj 
and I3] as the simplest model of intermediate statistics. 
It is defined in such a way that the probability of having 
N ordered levels Ei < E2 < ■ ■ ■ < En is 



PpiEi 



' En) ~ W \Ej+i - Ej 



(5) 



For this model all correlation functions are known when 
N —> 00 |2(]. In particular, the nearest-neighbor distribu- 
tion is 

P,{s)^A,s^e-^'^^'^%A,^^I±^, (6) 
the two-point correlation form factor takes the form 

(7) 



A>(t) = 1 + 2 Re — - 



(l-H27riT/(/3 4-l))'3+i-l ' 



2 



and the level compressibility which determines the asym- 
ptotic behaviour of the number variance is x = ^(0) = 
l/(/3+ 1) . The two-point correlation function is expressed 
through the Mittag-Leffler function For integer (3 = 
n — 1 it is a finite sum 



exp 



2m- 



(8) 



In Figs. Hand 121 a few examples of nearest- neighbor dis- 
tributions for quantum maps with different a and sym- 
metries are presented. For each matrix 100 different re- 
alizations of random phases were considered and in all 
figures the spectrum is unfolded to unit mean density. 
Prediction © with Q agrees very well for all a with 
numerical calculations. 




FIG. 1; The nearest-neighbor distribution for non-symmetric 
matrices with a = 1/3, 1/6, 1/9 and resp. = 202, 205, 206. 
Solid, dotted, and dashed lines are prediction © with resp. 
/3 = 2, 5, 8. Insert: the same but for M" where n is the 
denominator of a. For clarity pictures are shifted along the 
diagonal. Solid lines: the Poisson prediction p{s) = e~*. 

We sketch main steps leading to the above result. De- 
tails will be given elsewhere jlOj. First, for all a and 
N eigenphases of matrix Q are such that for all nearby 
pairs of eigenphases ipi and ip2 there is one and only one 
eigenphase (p such that ip + ^-na falls in-between and 
Lpi and there is one and only one ^p' such that the same 
is true for ip' — l-Ka 



< V + ^TTCk < </32 , '^i — 27ra < ip2 



(9) 



To check it let us consider instead of matrix Mkp © a 
new matrix 



Nkp = Mk. 



kp 



^27Ti{k-p-{-aN)/N 



(10) 



If Uk{i) and A(j) for j = \, . . . ^ N are eigenfunctions and 
eigenvalues of matrix Mkp {h.{j)uk{j) = J2p^kpUp{j)) 




FIG. 2: The nearest-neighbor distribution for symmetric 
matrices with a = 1/2,1/3,1/5,1/7 and resp. A'^ = 
201,202,201,202. For clarity pictures for a = 1/5 and 1/7 
are shifted horizontally by 1 and 2 units. Solid lines are pre- 
diction © with resp. /3 = 0,1/2,3/2,5/2. Insert: the same 
but for M". Solid lines: the super-Poisson distribution 11911 . 



eigenfunctions and eigenvalues of matrix Nkp- On the 
other hand, due to special form of matrix Mkp, matrix 
Nkp (|10|l is a rank-one perturbation of matrix Mkp 



kp 



M, 



1 - e 



2wiaN 



kp 



J'S'k 



N 



(11) 



then -^kU) = e^''^^'^Uk{j) and A'(j) 



A(j) are 



It is known (see e.g. [ll|) (and can be checked in this 
case) that for rank-one perturbations new eigenvalues are 
in-between the unperturbed ones. As both are known, 
the first inequality in follows. The second can be 
proved by similar arguments. 

These inequalities manifest the existence of long dis- 
tance correlations between eigenphases of matrix 
When condition ||2J) is fulfill, these long correlations in- 
duces a very particular short ordering. Namely, consider, 
instead of matrix M defined by Q its n}^ power, M", 
where n is the denominator of a. Of course, eigenvalues 
of this matrix are just the n^^ power of eigenvalues of 
matrix M but, as eigenphases are defined only modulo 
27r, new eigenphases are equal to Tnp[j) mod 27r. There- 
fore, even if between two eigenphases tpi and (p2 there 
were no eigenphases of matrix M, between rupi and rnp2 
in M" there will be, in general, a few new eigenphases 
coming from different numbers of rotations around the 
circle. Nevertheless, condition is sufficient to insure 
that in matrix in-between riLpi and np2 there exist 
exactly n — \ new eigenphases. 

To prove it we put all eigenphases of matrix M on the 
unit circle and divide it in n consecutive sectors of angle 
27ra = 27rm/n. Let Uk be the number of eigenphases in- 
side the fc'^ sector. Denote by Xk and yk the (positive) 
distances between angle 2Trm{k — l)/n (i.e. the begin- 
ning of the k^^ sector) and the two eigenphases closest 



3 



to this angle in such a way that Xk corresponds to an 
eigenphase in the fc**^ sector and yk to one in (A; — 1)*'' 
sector. Inequalities Q lead to the following inequality 
for all k El 



k+l 



yk){xk+i - Xk) < 0. 



(12) 



As in the (fc + 1)*'' sector it must exist eigenphases which 
after the shift by — 27ra come in-between rik eigenphases 
of the fc*'' sector, it follows that the number of eigen- 
phases in the (fc -f 1)'^ sector is at least nu — 1. The only 
possibility to get more eigenphases in the (fc-|- 1)*^ sector 
is connected with the positions of the last and the first 
eigenphases in the sector. Straightforward calculations 
plus (|12(l give recursive relations 

nfc+i = nfc - 1 + Q{xk+2 - Xk+i) + 6(a;fc - 2:^+1) . (13) 

Here Q{x) = 1 if a; > and 6(2;) = if x < 0. 

Now choose the beginning of the first sector at the 
position of an eigenphase (i.e. j/i = 0). From (|12|l it 
follows that X2 < xi and x„ < xi. Direct applications of 
()13|l show that for j = 2, . . . , n (with Xn+i = Xi) 



ij = rii + Q{xj+i — Xj) . 



(14) 



As ^fc — mA^ (because it corresponds to m full 

turns around the unit circle) one obtains 

n-l 

niN — nni + 1 + 0(xj+i — Xj) . (15) 

If mN = 1 mod n, all the 0-functions have to be equal 
to which leads to the following chain of inequalities 



{) < Xn < Xn-l < . . . < X2 < Xi 



(16) 



But xi is just the distance between a pair of two nearby 
eigenphases of matrix ^ and the true eigenphase corre- 
sponding to Xk is 2nm{k — l)/n + Xk- Distances of all 
other eigenphases from boundaries of sectors are bigger 
than xi. When matrix M" is considered, multiples of 
27r have to be ignored and one concludes that in-between 
two eigenphases of matrix M" coming from two nearby 
eigenphases of matrix M there are exactly n — 1 other 
eigenphases ordered as in (|16|l . When raN = — 1 mod 
n all the 8-functions have to be equal to 1 which leads 
to the same conclusion (but n — 1 eigenphases appear in 
the inverted order). In other words, when eigenphases 
of matrix M" are known and Q is satisfied, correct or- 
dering of eigenphases of matrix M corresponds simply 
to considering each v}^ eigenphase of M" (and dividing 
them by n). 

But the n"^ power of the classical map with rational 
a = m/n \s classically integrable (as the momentum is 
conserved) 



ja) 



Quantization of this map as in leads to a quantum map 
whose eigenphases are (po(fc) = —2nNjy^^iP{k/N + 
ja). If in the semiclassical limit N 00 classical be- 
haviour dominates, it is natural that local spectral statis- 
tics of A/" (for all A^) will be Poissonian as it is con- 
jectured for all generic integrable models 01 • 1^ seems 
that it is the case for non-symmetric matrices (with an 
arbitrary function f{k)) because all ipo{k) are different 
and small corrections from omitted terms are negligible. 
But for symmetric matrices F(—k) — F{k), consequently, 
ipoiN — k) = (/?o(fc) and these corrections are important 
as unperturbed eigenphases form degenerate pairs. Nev- 
ertheless, one can argue that each second eigenphase of 
(i.e. either odd or even in the ordered sequence) still 
form a Poisson sequence but degenerate levels are split. 
Let p{s) be the distribution of nearest-neighbor odd-even 
eigenphases with mean density 1. As odd or even eigen- 
phases of matrix Af " are assumed to obey the Poisson 
statistics (with density equals one half of the total den- 
sity) one gets a convolution equation 



whose solution is 



pis) 



p{l)p[s ~ l)dl 



-s/2 



(18) 



(19) 



mod 1 . (17) 



Higher correlation functions are obtained in the same 
manner and one concludes that the odd-even distribution 
is exactly the semi- Poisson distribution with /3 = —1/2 
(cf. ©). For brevity we refer to it as to the super- 
Poisson distribution because instead of level repulsion it 

shows level attraction (p(s) 00) [l^ . 

The above arguments demonstrate that for large A^ 
eigenphases of matrix Af" where n is the denominator of 
a have universal distribution (independent of a), namely, 
the Poisson distribution for non-symmetric matrices and 
the super-Poisson distribution (i.e. the semi-Poisson dis- 
tribution with [3 = —1/2) for symmetric ones (see inserts 
in Figs. HI and EJ. 

When A obeys Q eigenphases of M are obtained sim- 
ply by jumping on each n}^ eigenphase of matrix Af" 
which leads directly to l@J. In [l^ it was noted that for 
integer (3 the semi-Poisson distribution can be obtained 
from a Poisson sequence by considering each {(3-\- 1) level. 
It is interesting that this completely artificial mechanism 
appears naturally in the quantization of the map (0) 
(for symmetric matrices one jumps over not a Poisson 
sequence but a super-Poisson one). 

The described result applies only to local statistics in 
the universal limit of unit mean density and A^ — > 00. 
For long-range statistics (on distances of the order of 
A in this scale) the above-mentioned long-range corre- 
lations between eigenphases of matrix M become appar- 
ent. In Fig.|3|the two-point correlation function for non- 



4 



symmetric matrices with a = 1/20 and N = 801 is pre- 
sented. In the insert the correlation function till s = 5 is 
given together with the prediction © with n = 20. The 
very strong level repulsion (~ s^^) and prominent oscil- 
lations described by this formula agree very well with the 
numerical calculations. Nevertheless, for s close to aN 
one observes a big peak which is a clear signal of long- 
range correlations. At higher s new peaks (with smaller 
amplitudes) appear at multiples of this quantity. 
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FIG. 3: The two point correlation function for non-symmetric 
matrices with a — 1/20 and = 801. Insert: small-s be- 
haviour of this quantity. Solid line: Eq. with n = 20. 



We also investigate numerically eigenfunctions Uk{j) 
of matrix M. As expected for models with intermediate 
statistics they show fractal properties which manifest e.g. 
in a non-trivial power increase of participation ratios 



Rn 



(E 

k 



2qx 



N- 



where J^k — 1 (• ■ ■)j denotes an averaging 

over the spectrum. For example, for non-symmetric ma- 
trices with a = 1/2, I?2 ~ i's ~ -D4 w .4. 

In conclusion, we investigate unusual spectral statistics 
of random unitary matrices coming from a quantization 
of map 101. The most interesting case corresponds to 
a = m/n (when the map is an interval-exchange map) 
and mN = ±1 mod n [la|. When N ^ 00 local spec- 
tral statistics of these matrices tends to the semi-Poisson 
statistics with parameter /3 which can take arbitrary in- 
teger and half-integer values depending on the denomi- 
nator of a and the time-reversal invariance of the map. 
The level compressibility is non-zero (l^ . Eigenphases 
of M" for all N have the Poisson distribution for non- 
symmetric matrices and the super-Poisson distribution 
for symmetric ones. When mN = ±1 mod n eigenphases 



of M are obtained by jumping over each n^^ eigenphase 
of M" thus explicitly realizing the ad hoc mechanism 
proposed in Eigenfunctions have fractal properties 
in momentum representation. 

The model considered seems to be the first example of 
a new class of random matrix ensembles and its detailed 
analysis may lead to a better understanding of spec- 
tral statistics of pseudo-integrable plane polygonal bil- 
liards whose classical mechanics is described by interval- 
exchange maps. 

The authors are greatly indebted to J. Marklof for dis- 
cussing the paper jlj prior the publication and for numer- 
ous useful discussions and comments. One of the author 
(E. B) is thankful to O. Giraud for useful discussions and 
to Y. Fyodorov for pointing out Ref. Q . 
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